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Abstract 

Mutually conjugated synergetic schemes are assumed to address evolution of 
nonequilibrium self-organizing system. Within framework of the former, the sys- 
tem is parameterized by a conserving order parameter being a density, a conjugate 
field reducing to gradient of related flux, and control parameter, whose driven mag- 
nitude fixes stationary state. We show that so-introduced conjugate field and control 
parameter are relevant to entropy and internal energy, so that self-organization ef- 
fect is appeared as a negative temperature. Along the line of the conjugated scheme, 
roles of order parameter, conjugate field and control parameter are played with a 
flux of conserving value, and gradients of both chemical potential and temperature. 
With growth of the latter, relevant value of the entropy shows to decrease in su- 
percritical regime related to spontaneous flux-state. We proof that both approach 
stated on using density and conjugated flux as order parameters follow from uni- 
fied field theory related to the simplest choice of both Lagrangian and dissipative 
function. 
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1 Introduction 



Nowadays, complex phenomena type of nonequilibrium phase transformations, synergetics, 
self-organized criticality, etc. attract very much attention [1], [2]. That is stipulated not only 
applications to a vast variety of concrete problems, among which are pattern formation, mor- 
phgenesis, traffic flows at so one, but also due to lack of a general theory of nonequilibrium 
systems, in contrast to equilibrium statisical physics based on well-established concepts and 
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tools (see [3]). The latter are known to be stated on introduction of thermodynamic potentials 
and extremum principle [4]. It is quite clear that, within framework of the theory of nonequi- 
librium systems, these ingredients must keep in relations respected to equilibrium limit, along 
with Onsager-type relations between gradient of thermodynamic state parameters and con- 
jugated fluxes. This paper contains an attempt to build up such a scheme originated from 
synergetic concept. 

As is known, synergetics deals with far-from-equilibrium systems in a variety of fields, ranging 
from physics over chemistry and biology to sociology, where external driven affects arrive at a 
self-organization [5]. Formally, the latter is appeared as entropy decrease with the subsystem 
moving off equilibrium - in contrast to the second law of thermodynamics that is applicable 
to whole system, which keeps in equilibrium as rule. Being metastable in nature, such ordered 
state tends to a local minimum of synergetic potential, where it may be located during very 
long time intervals. Main difference of such states is to be further stationary, than equilibrium 
one. 

In general case, two types of stable states are possible: (i) with nonzeroth magnitude of an 
order parameter and (ii) with the same of a conjugated flux. It is appeared that the only first 
state is essential in the case when relevant order parameter is nonconserving. In this case, 
the synergetic approach addresses kinetic of nonequilibrium phase transition. However, in 
much more complicated case of a conserving order parameter a flux state could appear spon- 
taneously. In this connection, a problem arise to treat this state in analogy with synergetic 
theory of phase transition [6]. A method adopted in this work is stated on a generalization of 
such approach to a general theory of nonequilibrium systems based on axiomatic field theory. 

The paper is organized in the following manner. In Section 2, a synergetic scheme based on 
the famous Lorenz system is assumed to address evolution of self-organizing system param- 
eterized by a conserving order parameter, a conjugate field, reducing to gradient of related 
flux, and a control parameter, whose externally driven magnitude fixes stationary state of the 
system under consideration. We show, within adiabatic approximation, that in the case of 
nonequilibrium thermodynamic system so-introduced conjugate field and control parameter 
are relevant to the entropy and the internal energy, respectively. As a result, self-organization 
effect is expressed in appearance of a negative temperature, which magnitude decreases mono- 
tonically with growth of driven control parameter. Section 3 is devoted to study of flux state, 
whose parameterization is achieved by introducing, apart from a flux, of gradients of both 
chemical potential and temperature to play role of the conjugate field and the control param- 
eter, respectively. With growth of the latter, relevant value of the entropy shows to increase 
monotonically in subcritical regime and decrease in supercritical one related to spontaneous 
flux-state. Principle particularity of self-consistent equations used in Sections 2 and 3 is that 
these are postulated to be reduced to the form of familiar Lorenz scheme. The only justifi- 
cation of such a choice is that relevant steady state relations lead to the standard equalities 
of the equilibrium thermodynamics. Section 4 proofs that both approach stated on using 
density and conjugated flux as order parameters are followed from unified field theory for a 
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set of pseudo-vectors composed of: (i) density and flux; (ii) entropy and gradient of chemical 
potential; (iii) internal energy and temperature gradient. It is appeared that Lagrangian and 
dissipative function related to the Lorenz scheme are the simplest ones: the former contains 
quadratic and cubic terms only, whereas the latter - quadratic time derivative contributions, 
whose prolonging accounts dissipation processes. 



2 Density as order parameter 



We address a system with density p = p(r) and velocity V = V(r) distributed over space 
according to the following decompositions: 

p = Po [l + rj(r)}; V = V [1 + v(r)] , 1 = ^. (1) 

Here p , V = const are volume averaged values, rj(r) : v(r) are space-dependent dimensionless 
components. Relevant flux and energy 



Q^pV, £ = ^V 2 (2) 



(3) 



P 
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are determined by these components as follows: 

Q = Q [(l + 77)1 + q(r)], E = E [(l + V ) + 2q(r) + e(r)]; 

Qo = PoVo, E = 

2p 

Statistical state is fixed by internal components 

q=(l + ?7 )v, e=-^. (4) 

1 + T] 

A conserved value of density p is governed by the continuity equation 

t + VQ = 0, V-|, (5) 

where r, t are coordinate and time. With switching a dissipative sink — rj/r characterized 
with a relaxation time r, the continuity equation (5) takes the form 
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— + VbVq = , V = — . 6 

r p 

Introducing source term 



s = -AVq, A = V t, (7) 
one reduces this equation to the simplest form 

rrj = -7] + s, (8) 



where dot stands for time derivative. At stationary state (?) = 0) this equation shows coinci- 
dence of the source s and the density deviation rj. 

Let us find now an evolution equation for the source (7) related to the dimensionless flux q. 
This equation should be obviously reduced to Onsager relation 



q = -XVr] (9) 
at stationary state, where the source takes the form 

s = A 2 VV (10) 

On the other hand, here s = rj, so that the dependence rj(r) is governed by the Poisson 
equation 



\ 2 V 2 r] = r]. (11) 

To be subjected to above conditions, we postulate the evolution equation found in the fol- 
lowing form: 



t s s = -s + r]e. (12) 

Here, the first term in right-hand side is addressed relaxation process characterized with time 
t s , the second term is non-linear in the form due to a parameter e, being introduced to take 
magnitude e = 1 at stationary conditions. Time evolution of this parameter is supposed to 
be of relaxational type also, but if the autonomously evolving values r](t), s(t) go to zero with 
time t — > oo, the parameter e(t) is supposed to tend to a finite magnitude eo ^ 0, within an 
autonomous regime. We postulate the evolution equation for time-dependence e(t) of such a 
parameter as follows: 
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T e e = (e - e) - rjs, 



(13) 



where r e is relevant relaxation time. 

Being equivivalent to the Lorenz form (see Section 4), equations (8), (12), (13) obtained 
constitute the basis for self-consistent description of the evolving system with driven control 
parameter e and relaxation times t v , t s , r e . The distinguishing feature of these equations is 
that nonlinear terms, that enter Eqs. (12), (13), are of opposite signs, while Eq. (8) is linear. 
Physically, the latter means that on the early stage of the process under consideration, the 
dimensionless density deviation r\ is reduced to the source s defined Eq. (7). The negative 
sign of the last term in Eq. (13) can be regarded as a manifestation of Le Chatelier principle. 
Indeed, as we convince below, supercritical e-increase results in self-organization process, 
whereas the density deviation r\ and the source s in Eq. (13) tend to suppress the growth of 
this parameter. The positive feedback of r\ and e on s in Eq. (12) plays a fundamental part 
in the problem. As we shall see later, it is precisely the reason behind the self-organization 
that brings about the self-organization. 

In general case, the system (8), (12), (13) cannot be solved analytically, but in the simplest 
case, where r s , r e <C t v , the source s and the control parameter e can be eliminated by making 
use of the adiabatic approximation that implies neglecting of the left-hand sides of Eqs. (12), 
(13). As a result, the dependencies of e and s on r\ are given by 



Note that, under 77 is within the physically meaningful range between and 1, the parameter 
e is a monotonically decreasing function of the density deviation r/, whereas the source s 
increases with rj (at r\ > 1 we have ds/di] < and the self-organization process becomes un- 
stable). Such type behaviour permit to assume that the source s can be meant as the entropy 
increase caused by the density deviation during self-organization, whereas the parameter e is 
a related energy. Within the framework of such supposition, the dependence 




(15) 



following from Eqs. (14), arrives at the temperature 



T= — 

ds 



(16) 



as follows: 



So defined temperature increases monotonically with energy growth from magnitude T = 
at e = to infinity at the point e = €q/2. Here, the temperature T breaks abruptly to negative 
infinity and then increases monotonically again to initial magnitude T = at e = eo- This 
means that inside domain < e < e /2 the self-organization process is dissipative to behave 
in usual manner; contrary, within domain e /2 < e < e self-organization process evolves 
so, that an energy increase derives to entropy decrease, in accordance with negative value of 
temperature. 

Substitution of the first Eq. (14) into Eq. (8) yields the Landau-Khalatnikov equation 

" = ~W (18) 



with synergetic potential 

^ = ^ 2 -|ln(l + ^ 2 ). (19) 



For eo < 1, the ^-dependence of W is monotonically increasing and the only stationary value 
of i] equals zero, r\ e = 0. If externally driven energy eo exceeds the critical value, e c = 1, the 
synergetic potential assumes the minimum with nonzero steady state magnitudes of density 
deviation and entropy r] e = s e = V e o ~ 1 to be dependent on the driven energy eo and 
constant value of the stationary energy e e = 1. The temperature (17) takes the stationary 
magnitude 



Ve - 1 

T °-~ 

being negative within supercritical domain 
T decreases monotonically with the driven 
e = 1 to negative infinity at e = 2. This 
with growth of nonequilibrium power. 



(20) 

1 < eo < 2. Thus, the stationary temperature 
energy increase from the zeroth magnitude at 
means that self-organization stimulus increases 



3 Flux as order parameter 



The above consideration shows that the dissipative dynamics of the conserved system can be 
represented within the framework of the Lorenz model, where the density deviation rj plays a 
role of the order parameter, the source (7) is reduced to the entropy s being conjugated field 
to the order parameter and the internal energy e is the control parameter. The key point of 
this approach is that one postulates to address the density deviation rj(r,t) as fundamental 
field. However, a variety of the physical systems are known as well [1] - [3], where the flux 
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deviation q takes a fundamental part of the order parameter. To study relevant systems, we 
postulate the Lorenz system in the following form (cf. Eqs. (8), (12), (13)) 



r q q=-q + f, (21) 
r f f =-f + l f (qg), (22) 
r g g=(go-g)-l g (qf). (23) 



Here, r q , Tf, r g are relaxation times, If, l g are unit vectors along related directions and 
we introduce a constant vector g = golo, which physical meaning will be determined. At 
the stationary state, when q = 0, a force f reduces to the flux q. This means that in non- 
stationary case q^O this force takes the form of generalized Onsager equality (cf. Eq. (9)) 



f = -AV/x, (24) 

where a specific thermodynamic potential \x is introduced to be the chemical potential in 
usual case. Thus, we can conclude that the vector f is reduced to the thermodynamic force. 



To state a physical meaning of the third vector g, let us suppose, as above, that relaxation 
times are subjected to conditions Tf , r g <C r q . Then, within the adiabatic approximation one 
obtains: 



O (l()lq)lf5 „_ l + <?o(lglf)[[lglo]lq] q 2 ^ 



~ l + (l g l q )(l f l q )q 2 ' * l + (l g l q )(l f l q )q 2 ' 

where parentheses and square brackets denote scalar and vector productions, respectively. 
These equalities are reduced to form of the relations (14) if the vector g do not varies its 
direction, i. e., [l g lo] — 0. Moreover, it is enough for our aims to restrict ourselves addressing 
the simplest one-dimensional case. Here, the expressions (25) give the following state equation 
(cf. Eq. (15)) 



/ = VWo-</)- (26) 



Then, a derivative 

df 



(27) 



takes the form (cf. Eq. (17)): 
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(28) 



With g-growth within the domain < g < go, this derivative increases monotonically from 
— oo to oo taking the magnitude s — at g — go/2. This means physically that the quantity 
defined by equality (27) could be recognized as a specific entropy infinitely increasing with 
tending the control parameter go/ 2 < g < go to externally driven value go- Respectively, the 
control parameter g gets meaning of the temperature gradient taken with inverted sign: 

g = — AVT. (29) 



At the stationary condition, when q = 0, one obtains the steady-state magnitudes of flux, 
thermodynamic force and temperature gradient as follows: 

Qe = fe = \fg~o - 1, g e = 1- (30) 



On the other hand, steady-state specific entropy (cf. Eq. (20)) 



l -So/2 

S0 = 7^t (31) 

decreases monotonically from infinity to zero with the temperature gradient growth within 
supercritical domain 1 < g < 2. Thus, in accordance with above supposition, self-consistent 
evolution of flux q, thermodynamic force f and sign inverted temperature gradient g arrives 
at entropy decrease that means self-organization process at supercritical magnitudes of the 
control parameter g > 1. 



4 Field theory of self-organizing system 



According to our previous contribution [7], the Lorenz system is relevant microscopically 
to the bozon-fermion Hamiltonian of Dicke type. Phenomenologically, its consideration is 
possible only within supersymmetry field theory [8]. To take into account quite different 
commutation rules related to different freedom degrees, let us introduce some pseudo- vectors 




C = 




TC EE -CF, f EE -AV/i, 



(32) 
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to yield material fields 0, as well as mutually anticommuting components of conjugate fields 

— * — * 

J- and state parameters C (here, parameters 5, e have to be defined). The system behaviour 
is postulated to be governed by Lagrangian 



£ — £o + C-int 

that comprises of a bare part 



(33) 



Cq — (pJ-'Co, Co 



■ 

go 



and an interaction contribution 



(34) 



C-int — —<f)J~C. 

In analogy, dissipative function 



(35) 



n = n m + n f 

consists of a material component 



(36) 



Km = ^ (0)^0) 2 

and field gradient terms 



(37) 



n f = \{v,fif + l -(v c c)\ 

where the notions are introduced as follows: 



(38) 



- d - _ 
V z = r z — + 1; Tj, 



V Tf , 



(39) 



which are determined by prolonging derivatives V z to take into account Debay-type relaxation 
given by times f z with z = (f),C,J-. Then, Euler equations 



dC _ dC &R 
~dz~ ~ z d{V z z) ~ d{V z z) 



, Z = <j),J r ,C, 



(40) 
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arrive at basic system of differential equations for pseudo- vectors (32): 



Vff = C $-$C, (41) 
V c C = $f, 

— * 

where the second equality takes into account anticommutation relation for pseudo- vectors T 
and C. Equations (41) are reduced to the form 

rr] = -r} + (1 - e/e )s, 

T s s = -s + e r]-r)e, (42) 
r e e — —e + rjs; 



7- q q = -q + (i - s/g )f, 

Tf f = _f + l f (g q) - l f (5q), (43) 
r s 5 = -5 + l 5 (qf) 

that takes the standard appearance [5] if one takes into account determinations (32) and 
keeps the linear terms only in equations for rj and q. To reduce these systems to the form of 
equations (8), (12), (13) and (21) - (23), respectively, one needs moreover to account relations 

s = e -e, 5 = g -g; g = -AVT. (44) 



Formally, the systems (42), (43) differ from corresponding equations (8), (12), (13) and (21) 
- (23) to arrive at the synergetic potential 



instead of the dependence (19). At supercritical driven energy e > 1 this dependence takes 
the minimum at the point rjl = v /e^ — 1 related to stationary synergetic potential W(r] e ) = 
y'io — 1/2, whose value is lower than initial magnitude W(0) = €q/2. However, one keeps in 
mind that, within a phenomenological ideology, the expressions obtained have an asymptotic 
character to coincide only in vicinity of the critical magnitude e = 1. As a result, we may use 
the simplest of approaches suppressing nonlinear terms in the first equations of the systems 
(42), (43). 
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FIGURE CAPTIONS 

Fig. 1. The energy dependencies of the system temperatures: (a) nonstationary magnitude T 
versus ratio e/eo; (b) stationary temperature T versus eg. 

Fig. 2. The entropy dependencies on the temperature gradient: (a) nonstationary magnitude 
s versus ratio g/go; (b) stationary entropy so versus go. 
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